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Abstract 
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■ This talk describes the evolution of studies of chaos in Yane-Mills fields, 

(N 

gravity, and cosmology. The main subject is a BKL regime near the singularity 
t = and its survival in higher dimensions and in string theory. We also 
describe the recent progress in the search for particle-like solutions of the 



O 

m 



^ ■ Einstein- Yang-Mills system (monopoles and dyons), colored black holes and 

■ the problem of their stability. 

o 

o . 

O ■ I- INTRODUCTION 

"o 

5I1 I The great achievements of theoretical physics in the 20th century, Einstein's theory of 

general relativity and the Yang-Mills theory of non-Abelian gauge fields, have many common 

•rH . 

^ I properties: both are gauge theories and they are, in contrast to Maxwell's theory of the 
electromagnetic field, nonlinear. From this point of view, the chaoticity of the corresponding 
fields is not a surprise. On the other hand, there are many examples of stable solutions 
(solitons) of nonlinear field equations. Thus, the question of chaos in the equations of general 
relativity (GR) and the Yang-Mills (YM) equations is not a trivial and straightforward one. 
The evolution of our insight into chaos in the YM theory provides a good illustration for 
this nontriviality. 
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In the 1970s the problem of the integrabihty of the classical YM equations was actively 
studied. All attempts to find additional (to the trivial ones, such as energy) integrals of the 
YM equations in Minkowski space proved unsuccessful. This gave rise to a progran aimed 
at describing the YM dynamics not in terms of the potentials and fields, but rather in terms 
of the loop (string) variables In a Euclidean space, the famous classical solutions called 
instantons [@], were already known. But what was going on in (3+l)-dimensional space- 
time? Why are soliton-like solutions absent there? The answer is quite simple: The YM 
equations are chaotic even, as it is shown below, in the simplest conceivable limiting cases. 



II. CHAOS IN THE CLASSSICAL YM SYSTEM 

Consider a spatially homogeneous, sourceless YM field with potentials A°^(t){i, a = 1,2, 3) 
in the = gauge, for the gauge group SU(2). Such vector potentials can be considered 
as the relevant degrees of freedom in the infrared limit. They also may describe field con- 
figurations with a high gluon density (see ref. for details). With these approximations, 
the YM equatins are reduced to a discrete Hamiltonian dynamical system: 

^ - /(^ -4^' - = (1) 

with the conserved external and internal angular momenta: 

= e°'^^A\kl. (2) 

These quantities vanish for the sourceless field. Thus there exist, including the energy, seven 
integrals of motion, but the system is still not integrable. In the simplest nontrivial case, 
with X = gA\, and y = gA^ as the nonvanishing components, we obtain a system of two 
nonlinearly coupled oscillators with the Hamiltonian (the socalled x^y^ model ^): 

H=]-{pI+pI + xW) (3) 
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which exhibits strong chaotic behavior: most periodic orbits are unstable, one Lyapunov 
exponent is positive, etc. (see ref. for details). 

There are several mechanisms that can eliminate this chaoticity . One that is especially 
important here, is the coupling to a scalar field (Higgs field) with nonvanishing vacuum 
expectation value (0) = v. Then a term proportional to g'^v^{x'^ + y"^) is added to the 
Hamiltonian , which reduces the system to one of two weakly coupled harmonic oscillators 
when the vacuum expectation value v is large. However, at small enough v the chaos persists, 
and the transition between chaotic and regular behavior is controlled by the parameter 
g'^v'^/E where E is the energy density 

The chaoticity of the sourceless YM field theory was also established for the special 
case of spherically symmetric field configurations . The systematic stuy of the sourcesless 
classical lattice gauge theory by the Duke group also found evidence for spatio-temporal 
chaos in this system (see ref. for details). 

Concluding this brief survey of the study of chaos in the YM fields, I need to emphasize 
that these results demonstrate that the classical YM fields lack any special stable config- 
urations. All states are chaotic and no particular configuration dominates, in contrast to 
the Euclidean case, where the instantons make a dominant contribution to the functional 
integral. It is also noteworthy that the YM field is a dynamical system where quantum 
corrections are undoubtedly important. These corrections will lead to the elimination of 
most, but not all, properties that are characteristic of the chaotic YM system. This holds 
true for the YM quantum mechanics as well as for the quantum YM field theory ^j^. 

III. CHAOS IN GRAVITY AND COSMOLOGY 

A. BKL scenario of the approach to the singularity 

The study of chaos in general relativity (GR) has an even longer history than in the case 
of the YM field. Einstein himself believed that, due to the highly nonlinear self-interaction 
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of the gravitational field, the equations of motion of GR cannot be solved in general. This 
is the reason why he was so pleased when Karl Schwarzschild found the general solution 
with spherically symmetry. In fact, due to its different gauge group structure, gravity has 
more "capacity to resist" chaos than the YM field. From this point of view, and to exhibit 
the closest similarity to the spatially homogeneous YM field, let us consider the Einstein 
equations near the space-time singularity in the synchronous coordinate frame, where the 
metric g^i, can be expressed as a function of the synchronous time t only. 

Landau remarked (see ref. P]) that g = det^gf^u) tends to zero as t — 0, at some finite 
time interval, independent of the equation of state or the character of the gravitational field. 
The behavior of the metric ds"^ = dt^ — di"^ with {a,f3 = 1,2,3) di"^ = '-)apdx°'dx^ near the 
singularity at t = has the form 

de = t^P'dxl + t^P'dxl + t^P^dxl (4) 

with Pi being the eigenvalues of the matrix in the equation of motion: 

2 5 

Since A satisfies the trace relations tr(A) = tr(A^) = 1, the exponents obey the conditions 
X)i=iPi = YH=i'p'i = 1- Equation (||) is the so-called Kasner solution flOl corresponding to 
the Bianchi I geometry. Thus, in this simplest case, we have the regular flat, homogeneous, 
but anisotropic space with the total volume homogeneously approaching the singularity 

{dxidx2dx3 ~ if:). 

A generally chaotic oscillatory behavior was observed for the Bianchi IX geometry in the 



pioneering study of Belinski, Khalatnikov, and Lifshitz [^|Tl| of the Mixmaster universe [|12 
This study greatly stimulated the interest in the problem of chaos in GR. For the Bianchi 
IX geometry, a generalization of (H): 

de = {a^iaifj + h'^mamp + (?nanl3)dx°'dx^ (5) 

where £, m, n are unit vectors along the principal axes, the coefficients a, 6, c are functions of 
t. The evolution towards the singularity at t = proceeds via a series of successive oscilla- 



tions, during which the distances along two of the principal axes oscillate, while they shrink 
monotonically along the third axis (Kasner epochs). The volume V = IGn'^abc again drops 
approximately as t. A new "era" begins when the monotonically falling metric component 
begins to oscillate, while one of the previously oscillating directions begins to contract. This 
approach to the singularity, revealing itself as an infinite succession of alternating Kasner 
epochs, asymptotically takes on the character of a random process. 

In the study of the Mixmaster universe it is useful to employ its Hamiltonian formulation 
1^ , defined in the mini-superspace of the variables where fl is the volume and 



the P± are connected to the coefficients a ,b , c as follows: 

7„, = diag(a^6^c2)=e2^(e2^)^^. (6) 

Obviously, tr(/3) = 0. In this approach, the spatial scalar curvature V^(/?+,/5_) plays and 
essential role as the two-dimensional time-dependent potential in the Hamiltonian 

2H = -pI +pI+ p'_ + e-^^ {V{P^, p.) - 1) (7) 

with Pn,+,- being the momenta conjugate to the variables The time coordinate 

T is defined by dt = abcdr, where t is the synchronous time. For the Bianchi geometry, 

The evolution governed by epochs and eras here corresponds to bounces off the potential 
V which has three corners in the plane /?_) for fixed Q. The epochs correspond to 
bounces outward along the corner of the potential, while a new era begins with the change 
of a corner. Qualitatively, the picture resembles the behavior found in the x'^y'^ model 
of YM dynamics, where the motion along the channels formed by the equipotential lines 
(hyperbolae) corresponds to an "epoch" , and the change of the channel starts a new "era" 
(see ref. for details of the YM case). 

Concluding this subsection, we remark that a powerful tool for the study of chaos, the 
Lyapunov exponents A, sometimes does not apply in the case of gravity, leading to a con- 
tradiction between the analytical statements of the BKL type and numerical calculations. 



This is not surprising since A, measuring the rate of divergence between neighboring trajec- 
tories, depends on the choice of the time variable. E.g., A = in the BKL time variable 



T = J dt/abc, while A = oo in the metric time t One needs to have a covariant 

measure of chaos that does not depend on the choice of the time variable. For instance, 
the notion of integrability, i.e. the existence of first integrals of motion, is appropriate. 
Concerning the cases discussed above, it has been shown that the Bianchi type VIII and IX 
models do not possess such integrals [|T5[] . 



B. BKL problem with a scalar field 



It was shown [jl6| in 1973 (see also the recent paper [|T7|) that the addition of a massless 
scalar field eliminates chaos in the BKL problem. For illustration, consider the simpler 
case a scalar field conformally coupled to gravity in the Friedman- Robert son- Walker (FRW) 
universe. In the Hamiltonian version we have after proper rescaling of the scalar field 0: 

H{a, <p- Pa, p^) = ^ [~pI - ka" + + k4>^ + m^^^a^] (8) 

with a{t) the FRW scale parameter and k = ±1,0. (Compare (|^) with the Hamiltonian of 
the x'^y'^ model of the YM-Higgs system 0. The strange sign of the variables associated 
with the FRW scale factor a is specific to gravity.) Obviously, chaos does not appear for the 
massless scalar field (m = 0). For m ^ 0, one finds, in general, chaotic behavior due to the 
coupling between the "oscillators" a{t) and (f){t). 

If one adds to (^ the self-interaction of the scalar field (p and the contribution from the 
cosmological constant ^(A(/)^ + Aqo^) (where again is the rescaled scalar field and Aq is the 



rescaled cosmological constant [T^), then the Painleve analysis of the equation of motion 

d + ka — m'^acj)'^ — AqO^ = 

+ A;0 + m^a^^ + A0^ = (9) 

leads to the conclusion that the integrability test is passed only for a few special values of 
Aq, a, and k |jl8| (see also ref. [[I9i^ ). The case m = in (|) gives the reguar regime: 

6 



a(t) = 



2k 



(cosh Vkty 



''^{coshVkty^ 
A 



(10) 



C. Higher dimensions 

From the BKL scenario we have seen that, as one approaches the singularity at t = 0, 
spatial dimensions {d = 3) are continuously mixed by the dynamics in an infinite series 
of oscillations. This mixing, as well as the successive length of the oscillatory periods (i.e. 
the number of Kasner epochs in each era), have the character of a random process [jll]. 



All spatial dimensions stand on an equal footing. But what happens if one increases, in 
the spirit of the modern development of Kaluza-Klein models and string theory motivated 
models, the number of spatial dimensions? 

The problem can have importance in the context of the dimensional reduction that occurs 
with compactification. Consider the Kasner solution which generalizes the metric (^ for a 
space-time of dimension D = d + 1: 

i=l 

with Pi{x) time independent and satisfying Kasner 's conditions 

i:p^{x)=j:p.{xr = l, (12) 
1=1 1=1 

and time-independent one-forms uji{x) = tj{x)dx^ with det(£*) 7^ 0. 

Similar arguments as for d = 3 show that for ci > 10 the monotonic approach to t = 
according to (0) has the same degree of generality as the BKL oscillatory regime. They are 
based on the possible neglect of the spatial gradients (present in the spatial curvature) com- 
pared with the terms involving time derivatives, which are always of the relative magnitude 

in the vacuum Einstein equations, for any value of d. Thus, one requires: 

limt2Wi?^ = 0, (13) 



where is the dimensional spatial Ricci tensor. This leads to the condition for the 

"Kasner stability region": 



for i,j,k all different. When (|T^ contradicts the Kasner conditions ([T^), we unavoidably 



have the BKL scenario. As it was shown in ref. this indeed occurs for d < 10. However, 



for d > 10, the conditions ( p!^ are not contradictory and define an open region of the 
so-called "Kasner sphere" (|T2|). 

Thus the generalized Kasner solution ([Tl|) with monotonic power-law behavior of the 
spatial distances, becomes a general solution of the Einstein equations near a cosmological 
singularity when D = d + 1 > 11. Moreover, in this case the chaotic regime of BKL type 
becomes unstable and is replaced by the monotonic regime (|ll|) |2^-^ . Whereas for d < 10 



all spatial dimensions appear on an equal footing, as remarked earlier, the spatial directions 
are split into two groups for d > 10. One group is expanding, the other one contracting. 
Thus, some dimensions dominate over others as the singularity is approached, which may 
provide a "natural" explanation for the compactification to a lower number of dimensions 
near the singularity. 

D. Superstring theory and the BKL regime 

Superstring theories with their rich field content (graviton, dilaton, other scalars, vectors, 
p-forms, etc.) and large number of dimensions {D = 10, 11) may modify the BKL regime 



even at the classical level . In particular, the dilaton affects the formation of the Kasner 



sphere. Instead of ([l2|) , we have: 



d / d \ 2 d 

1=1 \i=i / i=i 



(15) 



Two additional sets of stability constraints, called electric and magnetic components and 
related to the contribution of the p-form fields Ap to the Einstein-dilaton equations appear. 
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in addition to the pure gravitational stability conditions ([T^). The positivity of these con- 
straints, together with ensures the stability of the generalized Kasner solution. The 
result of the study |25| is that for all superstring models there exists no open region of the 
Kasner sphere (|T5|) where all exponents are positive. In other words, the generic solution 
near t = for massless bosonic degrees of freedom exhibits the BKL behavior. 

Three ingredients were used to obtain this result: (i) inclusion of the p-forms in the field 
content of the theory; (ii) a sufficiently strong coupling of the dilaton to the graviton; and 
(iii) string dualities.^ The limitation of this study is the restriction of its arguments to the 
tree level. 

In conclusion, we remark that if the BKL-type regime survives in string theories, then the 
spatial inhomogeneity continuously increases toward a singularity p5[ ( "turbulent" universe 

IV. SOLITONS AND BLACK HOLES IN THE EINSTEIN- YM SYSTEM 

Recently, the interplay of gravity and the YM system was studied very actively. It is 
known that the sourceless YM equations have no particle-like solutions in Minkowski space 
29| , |30[| . The chaotic properties of YM dynamics ensure that no particular configuration 



dominates. In the pioneering paper numerical evidence was presented for the existence 



of a discrete family of globally regular solutions of the static Einstein- YM equations for the 
YM gauge group SU(2)). 

The existence of particle-like solutions of the Einstein- YM system is a result of the 
balance between the YM repulsive force and the gravitational attraction. The situation is 
similar to the addition of a scalar field to the free YM field resulting in the formation of the 



^We note that there exists an opposite statement in the literature |26| that the duality symmetry 
appears to be incompatible with the presence of the chaotic behavior near t — > 0. This statement 
is based on the qualitative Hamiltonian description of the Mixmaster models. 
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't Hooft - Polyakov monopole solution. Complementary to the study the existence of 
black hole solutions of the EYM system was established in ref. (colored black holes). 
This result contradicts the notion of the "baldness" of black holes (see e.g. |Q) and hints 
at the instability of the colored black hole solutions. 

Indeed, it was shown by various authors pi|-p^ that the particle-like solutions of the 
type as well as the non-Abelian SU(2) EYM black holes |^2|], are unstable for A > (flat 
space or de-Sitter space). Surprisingly, the situation is different for the anti-de-Sitter case, 
A < 0. It was shown that a continuous class of regular monopole and dyon solutions 
without cosmological horizon and "hairy" black hole solutions |^ exist in the EYM theory 
in any space that is asymptotically anti-de-Sitter (adS). For solutions without nodes in the 
non-Abelian magnetic field, the stability was established. The stability of nodeless dyon 
solutions is awaiting further study. 

The fact that for A < there exists a continuous set of such solutions, whereas for 



A > the set is discrete, invites the conjecture that the moduli space of solutions has a 
fractal structure. In conclusion, the existence of stable particle-like configurations may have 
important consequences for the evolution of the early universe if it ever passed through a 
adS phase. 
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